Ultrafast photo-doping and effective Fermi-Dirac distribution of the Dirac particles in Bi2Se3 
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We exploit time- and angle- resolved photoemission spectroscopy to determine the evolution of the out-of- 
equilibrium electronic structure of the topological insulator 612863. The response of the Fermi-Dirac distribu- 
tion to ultrashort IR laser pulses has been studied by modelling the dynamics of the hot electrons after optical 
excitation. We disentangle a large increase of the effective temperature (T*) from a shift of the chemical poten- 
tial (/i*), which is consequence of the ultrafast photodoping of the conduction band. The relaxation dynamics of 
T* and 11* are fc-independent and these two quantities uniquely define the evolution of the excited charge pop- 
ulation. We observe that the energy dependence of the non-equilibrium charge population is solely determined 
by the analytical form of the effective Fermi-Dirac distribution. 



The recent discovery of topological insulators (TIs) is re- 
newing the attention on the effects of spin orbit interactions 
(SOI) in solids, paving the road for the emergence of new 
quantum states of matter The SOI acquires particular 

relevance in the case of systems containing high-Z elements, 
leading to the lifting of the Kramers spin-degeneracy in bro- 
ken inversion symmetry systems, as described by Rashba Mr 
H], Dresselhaus (B\ and Rashba-Bychkov 13, [11]. There- 
fore, understanding the consequence of SOI is of primary im- 
portance also for future technological applications in spintron- 
ics. 

TIs are band insulators (semiconductors) where the conduc- 
tion and the valence band states have opposite parities and 
their energy ordering is inverted by the SOI 081 |l6j]. The most 
prominent feature of their electronic structure is the odd num- 
ber of spin polarized Dirac cones at the surface, connecting 
the opposite sides of the bulk band gap, resulting in topo- 




been discovered. Among these, Bi2Se3 represents a paradig- 
matic case, owing to the simplicity of its band structure char- 
acterized by a single Dirac cone 10 [gt] . 

In topological insulators, the spin helicity of the metallic 
surface state offers the unique possibility to support spin- 
current. Hence, spin-polarized charge distributions can be 
generated by circularly polarized light ll20[ |21|. The topo- 
logical protection of the linearly dispersing surface state is 
expected to strongly affect the scattering mechanisms of the 
Dirac particle, with respect to normal metallic states. In par- 
ticular, the different couplings to optical and acoustic phonon 
modes have been recently studied 122] and the optical exci- 
tation of long-lived electron population in the surface state 
might play an important role in forthcoming opto-spintronics 



equilibrium electronic properties of Bi2Se3, investigated by 
time- and angle- resolved photoemission spectroscopy (tr- 
ARPES). Although conventional ARPES, with its surface sen- 
sitivity, has proven to be effective and rich of information 
|[ll[3[la], the combined use of ultrashort laser pulses and angle- 
resolved photoemission experiments offers the unique possi- 
bility to explore the dynamical evolution of the charge carri- 
ers after an optical excitation, opening the door to the inves- 
tigation of mechanisms and states hidden to experiments in 
the frequency domain. After optical excitation, electrons ther- 
malize on a short timescale owing to the fast electron-electron 
interaction. The resulting hot electron gas is described by a 
Fermi-Dirac distribution with an effective electronic tempera- 



devices II22L 12311 . 



In this paper we report on the study of the out-of- 



ture T*{t) and an effective chemical potential iJL*{t) 11241 12511 . 
From the analysis of the temporal evolution of the Fermi- 
Dirac function at the Fermi wave- vector kp, we succeed in 
disentangling the large increase of T* from the ultrafast shift 
in the chemical potential, A/i*. 

Our experiment reveals that the dynamics of the non- 
equilibrium excited electron and hole populations strongly 
varies with the kinetic energy and the wave-vector {E, k). 
We develop a model based on the exponential relaxation of 
T*{t) and in the Fermi-Dirac function to quantitatively 
fit the measured dynamics. We find that the relaxation times of 
T*{t) and are /c-independent. The energy dependence 

of the charge population is then uniquely determined by the 
analytical form of the Fermi-Dirac distribution. Our model, 
owing to the universal properties of the Fermi-Dirac function, 
is more generally applicable to other tr- ARPES experiments 
on similar materials. 

High quality single crystals of Bi2Se3, in the form of 
platelets, were grown by the Bridgman technique. The Bi2Se3 
samples are n-doped due to atomic vacancies and excess se- 
lenium ]2]. The tr- ARPES experiments were performed at 
the T-ReX laboratory, Elettra (Trieste), with the use of a com- 
mercial Ti:Sa regenerative amplifier (Coherent RegA 9050), 
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FIG. 1: (color online) (a) ARPES band dispersion of Bi2Se3 at negative delay time along the TK high symmetry direction, (b) and (c) 
schematization of inter-band excitation channels between the conduction-conduction band (CB-CB) (b), and between the valence-conduction 
bands (VB-CB) (c). (d) and (e) Pump-probe tr-ARPES signal obtained as the difference between the images at positive delays (300 fs and 3.5 
ps) and negative delay (-500 fs). Red (blue) represents an increase (decrease) of the spectral weight, (f) Pump-probe difference EDCs at the 
Fermi wave-vector kp = 0.1 A~^, cut from panels (d) and (e). A shift of the chemical potential, Afi*, is experimentally resolved. 



producing 800 nm (1.55 eV) laser pulses at 250 kHz repeti- 
tion rate, with temporal duration of 60 fs. The samples were 
mounted on a He cryostat held at 100 K, and cleaved in situ 
in ultrahigh vacuum (2 • 10^^" mbar) at room temperature. 
The pump laser light was p-polarized, and the absorbed pump 
fluence was equal to 210 ixJ/cm?. The delay between pump 
and probe was introduced by modifying the optical path of the 
pump. Electrons were photoemitted by the i-polarized fourth 
harmonic at 6.2 eV, obtained by harmonic generation in phase - 
matched BBO crystals. The photo-electrons were analyzed 
and detected by a SPECS Phoibos 225 hemispherical spec- 
trometer, with energy and angular resolution set respectively 
equal to 10 meV and 0.3°, equivalent to 0.005 A^^ . The over- 
all temporal resolution was set to 300 fs, thus preserving the 
very high energy resolution required by the present experi- 
ment. 

Figure 1 (a) illustrates the band dispersion of Bi2Se3 at neg- 
ative delay time, i.e. before the arrival of the pump pulse (-500 
fs), along the VK high symmetry direction. The energy scale 
reports the measured kinetic energy, since the chemical poten- 
tial /i does not represent a good reference for the energy scale 
in the present experiment, as will be discussed in detail in the 
following. The linearly-dispersing topological surface state 
is clearly resolved. A black arrow points towards its apex, 
i.e. the Dirac point. The weaker signal within the Dirac cone 



at 1.7 - 1.84 eV corresponds to the bottom of the conduction 
band. These observations are consistent with previous ARPES 
results from conventional UV sources fllil. 

The intense pump pulses cause charge excitations between 
the occupied and unoccupied parts of the conduction band 
(Fig. 1 (b)), similarly to what observed in metals (CB-CB 
transitions) 11261 12711 . In the present case, the photon energy 
is larger than the band gap and this enables also inter-band 
excitations across the gap, from the fully occupied valence 
band (VB) to the partially unoccupied conduction band (CB) 
(VB-CB transitions) (Fig. 1 (c)). Figure 1 (d) and (e) show 
the pump-probe signal obtained as the difference between the 
ARPES images after and before the optical perturbation. Red 
(blue) represents an increase (decrease) of the spectral weight. 
At short delay time (h-300 fs) a large density of electrons pop- 
ulates the surface state and the conduction band above the 



chemical potential H28I1 . Such electrons result from intra- and 
inter- band transitions across the bandgap. However, the tem- 
poral resolution in this experiment is not sufficient to capture 
the formation of the nascent non-equilibrium electron distri- 
bution, which is expected to thermalize into a hot Fermi-Dirac 
distribution via electron-electron interaction, with a character- 
istic time shorter than 40 fs ll29ll . Figure 1 (e) illustrates the 
reduction in the pump-probe signal at larger delay times (h-3 .5 
ps). The density of excited charge carriers is reduced, as nor- 
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FIG. 2: (color on line) (a) EDCs at the Fermi wave- vector kp = 
0.1 at four different delay times (-700 fs, +300 fs, +3.5 ps and 
+7 ps). The fitting curves, in black lines, contain the FD distribu- 
tion, multiplied with a lorentzian which accounts for the surface state 
peak. The temporal evolution of the effective temperature T* and of 
the effective chemical potential /i*, as obtained from the the fit, are 
shown respectively in panel (b) and (c). A single decay exponential 
function well fits the temporal dynamics, with characteristic relax- 
ation time equal to tt ~ 2.5 ps and r,j ~ 2.7 ps, respectively. 



mally observed in metals 0261 12711 . 

Despite this similarity, we also observe a deviation from 
a purely metallic response. Such difference is well cap- 
tured by a close inspection in the pump-probe difference en- 
ergy distribution curves (dEDCs) at the Fermi wave-vector 
kp = 0.1 A^^ for various delay times fz^. Figure 1 (f) shows 
the two pump-probe dEDCs at +300 fs and +3.5 ps. We focus 
our attention on the zero-crossing point, which separates the 
charge-depletion and the charge-accumulation regions. The 
zero-crossing point clearly shifts more than 10 meV. We pro- 
pose that the modification of the zero-crossing energy can be 
explained by an ultrafast shift of the chemical potential, as 
recently reported in Bi2Se3 also by Wang et al. |22]. We in- 
terpret such an effect as the consequence of the photo-doping 
resulting from the inter-band excitations across the gap. In 
the case of a non-degenerate semiconductor (i.e. /i lying in 
the band gap) the optically excited electrons and holes must 
be considered as separated systems, each with a distinct ther- 
mal distribution, and with a different chemical potential de- 
pending on the charge density [30]. The thermalized electron 
and hole distributions relax towards equilibrium via radiative 
processes. Auger recombination or diffusion. A simpler relax- 
ation via diffusion was also proposed forp-doped Bi2Se3 123(1 . 
The electron-phonon scattering is unsuitable for recombining 
the excited electrons with holes, because the value of the gap 
is larger than the highest phonon energy (23 meV) ll23l IJH . 
Furthermore, for n-doped Bi2Se3, the band gap is acting as a 
bottleneck slowing down the relaxation of the photo-excited 
holes at high binding energies (in VB) and of the low energy 
holes (in CB). This results in a transient excess of charge car- 
riers in the conduction band, which is at the origin of the pro- 
posed ultrafast shift of the chemical potential. 



The emergence of the ultrafast chemical shift is further sup- 
ported by a quantitative analysis of the temporal evolution of 
the FD distribution function at kp. Figure 2 (a) displays the 
EDCs at some selected delay times (-700 fs, +300 fs, +3.5 
ps and +7 ps). The product of the FD distribution with a 
lorentzian, describing the surface state peak, is used as model 
function to fit the spectra |!28l]. The values of the effective 
temperature T* and of the effective chemical potential /i*, as 
obtained from the hot FD, are shown as a function of the delay 
time in Figure 2 (b) and (c) respectively. Their temporal evo- 
lution is fitted with a single decay exponential, with a proper 
rise-time (modelled by a step-function), convoluted with a 
gaussian to account for the temporal resolution. The char- 
acteristic relaxation times are tt ^ 2.5 ps and 2.7 ps. 
The T* relaxation is slightly slower than previously reported 
for p-doped Bi2Se3 (1.67 ps) [23]. This discrepancy might 
be ascribed to the different doping, but it could also be re- 
lated to the larger pump fluence used in the present case. The 
maximum value of A/i* is 12 meV, in agreement with the 
pump-probe dEDCs data. 

The creation of a long-lived population in the conduction 
band has been proposed to be the key-mechanism responsi- 
ble for the slow relaxation time in the topological insulators 
Bi2Se3 [23] and Bi2Te3 [3Z1 . In these models, the conduc- 
tion band acts as a charge reservoir for the topological surface 
state at lower binding energies. Such an effect was inferred 
from the difference in the temporal evolution of the tr-ARPES 
intensity at different positions in the band structure, I{E, k.t). 
Also in our investigation, the evolution of the surface state 
population at different binding energies along its linear disper- 
sion shows a peculiar energy dependent relaxation. However, 
a detailed analysis of the data suggests a different interpreta- 
tion. 

Figure 3 displays the I{E, k, t) signal, for six selected 
{E, k) positions below and above the chemical potential, as 
indicated in Figure 1 (e). We focus on the topological surface 
state and on the evolution of its out-of-equilibrium popula- 
tion. Similar trends are observed in the conduction band, as 
well. We average the intensity over a narrow energy window 
(less than 20 meV), and this provides us with a precise in- 
formation on the energy dependent temporal evolution of the 
excited electrons. This constitutes an important difference to 
a previous study on Bi2Tc3, in which large energy windows 
were selected, enclosing distinct bands [32]. The dynamics 
is fast for E » fjL (-6-), and apparently slows down when 
approaching fi (-4- and -3-). The behavior is opposite in sign 
but symmetric around /i (-1- and -2-). A similar observation 
was reported in previous studies of p-doped Bi2Se3, but the 
dynamics was analyzed only at a qualitative level [23]. In the 
following, a simple but effective model is proposed to account 
for this peculiar energy dependence of I{E, k, t). The spectra 
at the various kinetic energies are fit with the model function 
C{E, t), given by the following relation: 
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FIG. 3: (color online) Ultrafast evolution of the non-equilibrium charge population in the topological surface state at six energies, along its 
linear dispersion, as indicated in Figure 1 (e). The results of the fitting routine are displayed as black lines. The model function results from 
the intensity of the Fermi Dirac distribution at that energy at an effective temperature T* and with an effective chemical potential which 
relax exponentially in time as shown respectively in Figure 2 (b) and (c). 
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At time t the intensity (^{E, t) is the result of the convolu- 
tion of a gaussian G, describing the temporal resolution, with 
a Fermi-Dirac distribution. The latter is defined by a time- 
dependent effective temperature T* (t) and an effective chem- 
ical potential which relax in time as discussed previ- 
ously. O is a step function which accounts for the rise-time. 
The intensity C has a weak dependence on the kinetic energy 
E, which we attributed to matrix element effect. C accounts 
for the variation in sign in the signal between accumulation 
and depletion, when moving from above to below the chemi- 
cal potential. Finally, (f) represents a background. In our model 
the relaxation time is uniquely determined by tt and r^, and 
the apparent energy dependence is conveyed by the Fermi- 
Dirac distribution. The ({E, t) function is used to fit all the 
spectra (from -1- to -6-). The results of the fit are reported as 
black lines in Figure 3, showing a remarkably good agreement 
with the experimental data. 

A microscopic description of the temporal evolution of 
the charge population is complex, since several de-excitation 
mechanisms must be considered, especially for the recombi- 
nation of the photo-carriers across the band gap. Previous 
tr-ARPES study on p-doped Bi2Se3 and Bi2Te3 Q 
analyzed in detail the time evolution of the tr-ARPES inten- 
sity I{E, k, t). In these studies a delay between the dynam- 
ics of the surface state and the conduction band charge car- 
riers was revealed. This delay was interpreted as the result 
of a multi-steps relaxation process of the electrons excited in 
upper branches of the conduction band. Therefore, the con- 
duction band was proposed to act as a charge reservoir for 



a long-lived electron population in the Dirac cone. The Fig- 
ure 3 data analysis proves that the long lasting electronic dy- 
namics in the photodoped Bi2Se3 surface state is governed 
by the evolution of the Fermi-Dirac distribution. Recently, 
also Wang et al. focused their attention on the time evo- 
lution of the out-of-equilibrium Fermi-Dirac distribution, in 
Bi2Se3 12211 ■ The authors, by studying as a function of the 
temperature the response of Bi2Sc3 with varied stoic hiome- 
tries, were able to disentangle the role played by optical and 
acoustic phonons in the scattering mechanisms between the 
conduction band and the surface state. From the evolution of 
T* and /i* they pointed out that at low temperature the cou- 
pling between the bulk and the surface state is suppressed. 
However, in the present experiment, because of the tempera- 
ture (100 K), it was not possible to separate the acoustic and 
optical phonon cooling. Nevertheless, the scope of our work 
is to interpret the temporal evolution of the tr-ARPES inten- 
sity I{E, k,t). This is properly modelled by a time-dependent 
effective Fermi-Dirac function, uniquely determined by the 
effective electronic temperature T* and the effective chem- 
ical potential /i*. Moreover, our findings indicate that this 
pump fluence, and the consequent large increase of the ef- 
fective electronic temperature, might hinder the formation of 
light-induced spin currents at the Dirac cone. This suggests 
that the low fluence regime must be investigated in detail in 
order to reveal some of the non-trivial aspects of the carrier 
dynamics. 

In summary, we investigated in detail the time-dependent 
electronic thermal distribution after optical excitation in 
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Bi2Se3. We find that, at all the delay times, the electrons 
in the surface states are thermalized. Their distribution can 
be approximated by an effective Fermi-Dirac function, whose 
temperature T* and chemical potential fi* relax in time with 
a decaying exponential behavior. The fast electronic thermal- 
ization, below the experimental resolution (300 fs), is ensured 
by the characteristic time of the electron-electron interaction 
(< 40 fs 12% '). The variation in the effective chemical po- 
tential is interpreted as the result of the ultrafast photodoping 
of the conduction band, similarly to what recently reported by 
Wang et al. [22]. The relaxation time of T* and fi* are respec- 
tively equal to tt ~ 2.5 ps and ~ 2.7 ps. The former is 
compatible with a mechanism of energy relaxation to the lat- 
tice mediated by the electron-phonon interaction. The latter is 
related to the relaxation of the excess charge in the conduction 
band via diffusion. 

We propose a model which relies only on the inherent prop- 
erties of the Fermi-Dirac distribution function, therefore it 
provides a viable interpretation for similar energy dependence 
of the electronic dynamics in a different material |32l| or in 
p-doped Bi2Se3 ll23ll . The effect of the topological protec- 
tion during the relaxation of the Dirac cone charge population 
may provide the clue for designing new opto-spintronic de- 
vices lEoll . Our model paves the way for future investigations 
between the response of topologicaUy protected and conven- 
tional semiconductors. 
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